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We present an extension of the HTL resummation tech¬ 
nique to non-equilibrium situations, starting from the real 
time formalism in the Keldysh representation. As an exam¬ 
ple we calculate the HTL photon self energy, from which we 
derive the resummed photon propagator. This propagator is 
applied to the computation of the interaction rate of a hard 
electron in a non-equilibrium QED plasma. In particular we 
show the absence of of pinch singularities in these quantities. 
Finally we show that the Ward identities relating the HTL 
electron self energy to the HTL electron-photon vertex hold 
as in equilibrium. 


I. INTRODUCTION 


Here we use the following notation: K = (fcg, k), k = |k|. 
In equilibrium the distribution function fg is given by 
the Bose distribution nsiko) = l/[exp(|fco|/T) — 1]. In 
non-equilibrium the equilibrium distributions should be 
replaced by Wigner functions Jb = fB{\ko\,k,x) 0. 

A particular useful representation of the propagators in 
the RTF is the Keldysh representation, where retarded, 
advanced, and symmetric propagators are defined as lin¬ 
ear combinations of the four components of (||), of which 
only three are independent ||^: 

Dr = Dll — Di2, 

Da = Dll — D 21 , ( 2 ) 

Dp = Dll + D22- 


Perturbative QCD at finite temperature and density 
can be used to describe properties of an equilibrated 
quark-gluon plasma (QGP). However, using only bare 
propagators and vertices infrared divergent and gauge 
dependent result are obtained. These problems can be 
partly avoided by usin^he hard thermal loop (HTL) re¬ 
summation technique (y. In this way consistent results, 
which are complete to leading order in the coupling con¬ 
stant and gauge independent, are found. At the same 
time the infrared behavior is improved by the presence 
of effective masses in the HTL propagators. This method 
has been applied to various quantities of the QGP, such as 
damping rates of quarks and gluons, energy loss of ener¬ 
getic partons, thermalization times, viscosity, and photon 
and dilepton production [||. 

The HTL method has been derived within the imag¬ 
inary time formalism, which is restricted to equilibrium 
situations. However, in relativistic heavy ion collisions 
the formation of a pre-equilibrium parton gas has to 
be expected. Indeed it is not clear whether a complete 
thermal and chemical equilibrium will be reached in the 
evolution of the fireball at all ||]. Therefore it is desir¬ 
able to extend the HTL resummation technique to non¬ 
equilibrium. For this purpose we have to start from the 
real time formalism (RTF) [Q. Similar investigations 
have been performed by Baier et al. [^. 

In the RTF formalism propagators are given by 2 x 2- 
matrices . For example the scalar propagator reads 


DiK) = 


0 


0 

-1 


2TTiS{K^ 


( Ib 0{-ko) + fB\ 

V Oiko) + fB fB )- 


m^) 


( 1 ) 


For example, the bare scalar propagator in the Keldysh 
representation is given by 

Dr{K) = —^ 

+ i sgn(fco)e 

^ K2-m2-zsgn(A:o)e’ 

Df{K) = -2TTi (1 -f 2 /b) 6{K‘^ - rri^). 

The distribution functions only appear in the symmetric 
component Dp, which is of particular advantage for the 
HTL diagrams where the terms containing distribution 
functions dominate. 


II. HTL PHOTON SELF ENERGY 

Before we are going to calculate HTL’s in non¬ 
equilibrium we would like to note here that we do not 
aim to describe the thermalization of a partonic fireball 
starting far from equilibrium. As discussed in the liter¬ 
ature (see e.g. [^) this problem requires the use of non- 
perturbative methods beyond the HTL perturbation the¬ 
ory. We have situations in mind where we can use qua¬ 
sistatic distributions, i.e., where the equilibration is slow 
compared to the process under consideration. Such a sit¬ 
uation might be the case for the chemical equilibration of 
the parton gas in relativistic heavy ion collisions after a 
thermal equilibrium has been achieved rapidly 0. Also 
anisotropic momentum distributions taking account of 
the longitudinal and transverse expansion might be con¬ 
sidered 1^. At least as long as we are close to equilib¬ 
rium the use of quasistatic distributions within the HTL 


1 






method, based on the weak coupling limit g ^ 1, is jus¬ 
tified by the fact that the equilibration time r ~ 
is large compared to the HTL time scale {gT)~^. 

We have chosen a non-equilibrium QED plasma as an 
example for the use of the HTL method as the computa¬ 
tion of the HTL photon self energy is much easier than 
the one of the HTL gluon self energy. After all, the re¬ 
sults for these self energies differ only by a constant color 
and flavor factor. 

The HTL photon self energy is given by the one-loop 
diagram shown in Fig.l, where the internal momenta are 
hard compared to the external. In equilibrium this can 
be achieved in the weak coupling limit by considering 
internal momenta of the order of the temperature T or 
larger and external of the order gT or smaller. In non¬ 
equilibrium we might replace the temperature by the av¬ 
erage momentum following from the distribution function 
under consideration. 


K 



Here we restrict ourselves to the longitudinal compo¬ 
nent of the self energy, i.e. Hl = Hoo- Using the Keldysh 
representation and S{K) = 1^A{K) for the fermion prop¬ 
agators, we find 


U^{P) = -2ie^ 


d^K 


(qofco + q • k) 


Af{K - P)An{K) + Aa{K - P)Ap{K) 


(4) 


In the HTL approximation (AT ^ P) the retarded (Ha; = 
Hii -I- ni 2 ) and advanced (Ha = Hu -1- n 2 i) self energies 
are given then by 


where the constant A is given by 

^ dk P fpjk, x) [I - fpjk, x)] 

dkk fF{k,x) 

In equilibrium this constant reduces to A = 2T. 


III. RESUMMED PHOTON PROPAGATOR 

The resummed photon propagator is constructed from 
the Schwinger-Dyson equation shown in Fig.2. 


p 

= /\/\/\ + 


FIG. 2. Resummed photon propagator 



The Schwinger-Dyson equation for the resummed re¬ 
tarded and advanced propagators for longitudinal pho¬ 
tons reads 


^*R,A — ^R,A + Df,a^r,aD*r,a 


leading to 

D*iAip) = 


. 2 ^ Po,Po+P±ie 


p + I - — In 


2p po-p±ie 


(9) 


( 10 ) 


where the HTL photon self energy has been used. Apart 
from m-y, defined by (H), this expression is identical to 
the equilibrium one. 

The Schwinger-Dyson equation for the resummed sym¬ 
metric propagator is given by 

D*f = Dj:. + + D^U^D*a + Dj^U^D*A. 

( 11 ) 


This equation is solved by the following ansatz: 


K,AiP) = 



Po+p±ie \ 

po-p±ieJ ’ 


(5) 


where the effective photon mass is given by an integral 
over the non-equilibrium Fermi distribution fp 

dkk fpik^x). (6) 

In equilibrium this expression reduces to = e^T^/9. 

The symmetric self energy (Hi? = Hu -|- n 22 ) in the 
HTL limit can written as 



n^(P) = 2iA 


I^^r{P) 

Po 


(7) 


D*Up) = (1 + 2/b) sgn(po) [D*Up) - D*a{P)] 
+{n^(p) - (1 + 2 /b) sgn(po) iKiP) - ni(P)]} 
xD*^^{P)D*^{P). (12) 

The product of a retarded and an advanced propagator 
at the same momentum P may lead to a pinch singular¬ 
ity [^. In equilibrium, however, the term in the curly 
brackets vanishes as a consequence of the Kubo-Martin- 
Schwinger boundary condition or the principle of detailed 
balance. Then the full symmetric propagator is given by 
the first line of (^^, which is free of possible pinch prob¬ 
lems and reflects the dissipation-fluctuation theorem . 

The following relations can easily be shown to hold in 
general: 
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D*Up)D*a{P) 


D*l^P)_D*l^P) 
2iImIl^{P) ’ 
2iImIl^{P). 


(13) 


Together with the expression for the symmetric self 
energy we find from these relations that the symmetric 
non-equilibrium HTL resummed propagator can be writ¬ 
ten as 


D*Up) = - [D*Up) - D*^AiP)]- (14) 

Po 


Hence the symmetric propagator does not suffer from a 
pinch singularity even in non-equilibrium. This obser¬ 
vation also holds for p^—p^ > 0, where ImWR ~ e as 
following from (|^), because of the cancellation of ImH/j 
according to (|), (|T|, and H). 


IV. ELECTRON INTERACTION RATE 

As an example for the application of the non¬ 
equilibrium HTL resummation technique we consider the 
interaction or damping rate of an energetic electron in a 
non-equilibrium QED plasma. This quantity is one of the 
mostly studied within the HTL method. (For references 
see Q].) The HTL method does not lead to an infrared 
finite result in this case since there is not enough screen¬ 
ing in the transverse part of the HTL photon propagator. 
After all we want to study this example, since we are only 
interested in the comparison between the equilibrium and 
the non-equilibrium case. 

The interaction rate is defined by jl^ 

r(p) = -TV “ fpiP.x)] tr[f’ ImY,R{po =p,p)], (15) 

2p 

where 12r is the retarded electron self energy. Using only 
bare propagators the lowest order contribution to ImYiR 
comes from the two-loop diagram of Fig.3. 



ergy to the electron interaction rate within naive perturbation 
theory 

This diagram contains a dangerous pinch term coming 
from the product of a retarded and an advanced bare 
propagator at the same momentum Q |I4[ : 

Dr(Q)Da{Q) = 

- [sm?- (16) 


In equilibrium this term vanishes due to detailed balance. 

Applying the HTL method to the electron interac¬ 
tion rate results in using the one-loop diagram of Fig.4 
containing a resummed photon propagator for the elec¬ 
tron self energy |^J^. Owing to the non-zero imagi¬ 
nary part of the HTL photon self energy appearing in 
the resummed photon propagator this diagram exhibits 
an imaginary part. Since there is no pinch singularity 
in the HTL resummed photon propagator this diagram 
is also free of pinch singularities. In other words, the 
use of the HTL method to lowest order reduces two-loop 
diagrams, containing possible pinch singularities, to one- 
loop diagrams with resummed propagators, where pinch 
terms are absent. 



FIG. 4. Lowest order contribution of the electron self en¬ 
ergy to the electron interaction rate within the HTL method 


Following the equilibrium calculation [p|jl^] , we end up 
with the following result for the non-equilibrium interac¬ 
tion rate of an energetic electron: 


neq 


ip) ^ 


e^A , const 


Att 


In 


(17) 


The logarithm in this expression reflects the presence of 
an infrared singularity, which has been regularized by 
assuming an infrared cutoff of the order e^T. The con¬ 
stant under the logarithm cannot be calculated within 
the HTL method. As expected this result reduces to the 
equilibrium one |^, if we replace A by 2T.. 

Pinch terms may arise, however, in higher order contri¬ 
butions of the HTL perturbation theory, such as the one 
shown in Fig.5, where the product of a retarded HTL and 
an advanced HTL propagator at the same momentum Q 
shows up. 

Resumming the photon self energy H, which contains 
HTL electron propagators and HTL vertices, leads to a 
dressed propagator beyond the HTL resummation. The 
symmetric component of this propagator reads 


D**UQ) = -[D^*^R{Q)-D**^Am ( 18 ) 

Qo 

+{n^(Q) - - [n|i(Q) - n^{Q)]} d**UQ) d**a{Q)- 

qo 

If Hp’ is proportional to ImflR as in (|^), this propa¬ 
gator would not contain a pinch term. 
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FIG. 5. Photon self energy n 


V. WARD IDENTITIES 

Finally we will show that the HTL electron self en¬ 
ergy and the HTL electron-photon vertex fulfill the same 
Ward identities in non-equilibrium as in equilibrium. For 
this purpose we write the one-particle-irreducible vertex 
functions as |p^ , |T5| 

Tfi = ^112 + r2ii -I- r2i2 

— ri21 — Fi 22 — r221 — F222), 

Tfli = + rii2 -I- ri2i -I- F 122 

—r2ii — F212 — r22i — F222), 
r_Ro = 2^^111 + ri2i -I- r2ii - 1 - r22i 

— rii2 — F122 — r2i2 — F222)- ( 19 ) 

Here only the three retarded combinations of the seven 
independent components of the vertex in the Keldysh 
representation are shown. Using the HTL approximation 
for the vertices and the electron self energy, we can show 
the following Ward identities (without performing the 
loop integrations explicitly) 0: 

A^F^(-Q, P, K) = -fe(S^(P) - Sh(Q)), 

= -fe(S^(P) - S^(Q)), 

K^^Roi-Q^ P^ K) = -*e(Efl(P) - Efl(Q)), (20) 
which hold in equilibrium as well as non-equilibrium. 


VI. CONCLUSIONS 

In the present work we have extended the HTL re¬ 
summation technique to non-equilibrium situations as¬ 
suming quasistatic distribution functions. We are not 
able to treat the evolution of the system towards equi¬ 
librium, but can describe properties of a relativistic non¬ 
equilibrium plasma, where the equilibration is slow com¬ 
pared to the process under investigation. Such a situa¬ 
tion might be realized for example in a chemically non- 
equilibrated QPG produced in the fireball of a relativistic 
heavy ion collision. 

Starting from the RTF we have demonstrated the 
usefulness of the Keldysh representation for the HTL 
method in equilibrium as well as non-equilibrium. 


The non-equilibrium HTL resummation technique dif¬ 
fers from the equilibrium one only by two aspects: 

1. There are new effective masses defined by integrals 
over the non-equilibrium distributions appearing in the 
HTL self energies. 

2. There is a constant factor A, also defined by inte¬ 
grals over the distributions, in the symmetric self energy 
and the resummed symmetric propagator, which reduces 
to 2T in equilibrium. 

We have shown that there are no pinch terms in the 
HTL resummed propagators. Therefore pinch singular¬ 
ities are absent in quantities calculated from one-loop 
HTL self energies, such as the electron interaction rate. 
A similar result was found in the case of the photon pro¬ 
duction rate j|] and the collisional energy loss of energetic 
quarks |]l7|| . 

Beyond leading order in the HTL perturbation theory 
we expect also the absence of pinch terms using a higher 
order resummation scheme. 

Finally we have given the Ward identities for the HTL 
electron self energy and the HTL electron-photon vertex, 
which hold in equilibrium as well as non-equilibrium. 
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